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Abstract 

In an Endex reactor endothermic and exothermic reactions are directly ther- 
mally coupled and kinetically matched to achieve intrinsic thermal stability, 
efficient conversion, autothermal operation, and minimal heat losses. Ap- 
plied to the problem of in-line carbon dioxide separation from flue gas, Endex 
principles hold out the promise of effecting a COvcapture technology of un- 
precedented economic viability. In this work we describe an Endex Calcium 
Looping reactor, in which heat released by chemisorption of carbon dioxide 
onto calcium oxide is used directly to drive the reverse reaction, yielding a 
pure stream of CO2 for compression and geosequestration. In this initial 
study we model the proposed reactor as a continuous-flow dynamical system 
in the well-stirred limit, compute the steady states and analyse their stability 
properties over the operating parameter space, flag potential design and op- 
erational challenges, and suggest an optimum regime for effective operation. 
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1. Introduction 



In this paper we introduce a method of reactive endothermic-exothermic 
coupling, known as an Endex system [l|, for the in-line separation of car- 
bon dioxide from flue and fuel gas emissions. Endex principles are applied to 
obtain a highly efficient modification of the Calcium Looping separation tech- 
nique, in which calcium oxide (CaO, or lime) is used to scrub CO2 from a flue 
gas [3]. The Endex process involves direct thermal coupling of a carboniser 
reactor segment, in which the lime sorbent reacts with CO2 to produce cal- 
cium carbonate (CaCC^) and a scrubbed gas effluent stream, and a calciner 
reactor segment to which the loaded sorbent is transported and where the 
sorbent is regenerated and reinjected into the carboniser segment, with pro- 
duction of a pure CO2 gas stream. The same methodology can be applied 
to other sorbents and to fuel gas mixtures such as syngas. Since the Endex 
Calcium Looping reactor is a thermoreactive system an important first step 
in a project to build and operate a demonstration plant is to assess the ther- 
mal stability of the system. In this work we model the reactor as a coupled 
dynamical system in the well-stirred limit and map the linear stability of the 
steady states over the operational parameter space. 

Technologies for scrubbing carbon dioxide from flue gas to reduce green- 
house gas emissions must satisfy the following criteria, at least, in order to be 
accepted by society and implemented by industry: they must (1) effect C0 2 
capture for less than ~US$30 per tonne with an increased cost of electricity 
of less than about 15%, (2) preferably be retrofittable to existing fossil-fueled 
plants, (3) not cause additional environmental harm, and (4) operate safely. 
The best available technologies may not be able to comply with all of these 
requirements jij. There is a drive to develop second generation technologies 
which have been purpose-designed to meet the stringent demands of CO2 
capture, rather than adapt existing technologies. 

Calcium Looping was first proposed by Heesink and Temmink in 1994 
ji[ for post-combustion C0 2 emissions reduction, although its use on an 
industrial scale for removal of C0 2 from syngases dates from the 1960s [sj. 
Shimizu et al. [3] proposed that oxyfuel combustion of additional fossil fuel 
be used to drive the high temperature calcination process. Abanades et 
al. jsj showed that the C02-lime chemical bond strength is such that the 
thermal energy required is about 30-40% of the thermal energy of the plant 
that produces the flue gas. This additional consumption of energy results in 
higher CO2 production, but this is offset by use of the heat liberated from 
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the carboniser at about 873-1023 K to produce additional power. 

The major drawback of Calcium Looping is the rapid loss of carbonation 
capacity due to sintering of the lime sorbent. Sintering occurs at the high 
temperatures used in the conventional calcination cycle, about 1123-1223 K 
in an atmosphere of carbon dioxide 0, [lo[ and, more importantly, reactive 
sintering results in irreversible loss of the mesopore surface area due to pore 
filling [nj. 

Sintering was recognised as a problem by Curran et al. 0] and more 
recently by Silaban and Harrison [l2| and has been the subject of intensive 



research since. For example, Abanades and Alvarez [13] showed that the 
sorption capacity degrades from a capacity of about 80% of the theoretical 
limit of 0.78 kg of C0 2 per kg of CaO to about 17% after about 10 cycles, 



and reactivation studies were carried out by Fennell et al. [14| and Manovic 



and Anthony [151 ] . Even though MacKenzie et al. [16j showed that CO2 
capture using Calcium Looping with sintering meets criteria (l)-(3) above, 
sorbent sintering remains the largest barrier to adoption of Calcium Looping. 
In the Endex Calcium Looping process described in the present paper the 
sintering problem is neatly sidestepped by inversion of the usual temperature 
difference between carboniser and calciner and by using light carbonation 
over fast cycling. 

The Endex approach to thermal control and heat recovery was first de- 
scribed and analysed in mathematical terms by Gray and Ball 0, 0j- It 
involves using the heat generated by an exothermic reaction to drive an 
endothermic reaction directly, in real time. The reactions are thermally cou- 
pled via direct heat transfer or mass transfer or both, and for the reactions 
to effectively "see" each other they must be matched kinetically Since the 
kinetic parameters (activation energies and pre-exponential factors) of a se- 
lected Endex couple do not match in general, kinetic matching is achieved (to 
an acceptable approximation) by manipulation of residence times. Several 
studies of specific Endex systems have been carried out since with promis- 
ing results, although none has referred to the original works of Gray and 
Ball [l|, Q • Terms such as "recuperative coupling" and "co-current thermally 
coupled reactor" are used for Endex systems and principles in those works, 
which are reviewed in a recent study of a thermally coupled reactor [17J . 



Sceats [18| and Sweeney and Sceats 19] proposed that an Endex con- 
figuration could be used in which the temperature of the calciner is held 
below that of the carboniser by thermally coupling the reactor segments. 
This configuration is counter-intuitive because it is usually assumed that 
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the endothermic, bond-breaking, process in the calciner would occur at a 
higher temperature than that of the carboniser, and all previous approaches 
to Calcium Looping are based on this conventional configuration. Inversion 
of the temperature difference between carboniser and calciner is achieved by 
control of the pressure. Operation using this new configuration has three 
consequences: the lower temperature of the calciner reduces thermal sinter- 
ing, partial carbonation reduces the loss of surface area by pore-filling, and 
the unreacted sorbent promotes heat transfer. The intent was to greatly 
reduce the impact of thermal and reactive sintering, and remove this bar- 
rier to adoption of Calcium Looping. The most desirable property of the 
Endex configuration is that CO2 separation can be achieved, in principle, 
without the need for additional heat. Sweeney and Sceats [19| conclude that 
the Endex configuration has the potential to significantly reduce the cost of 
CO2 capture below that of current technologies, principally because of this 
property. 

A schematic of the Endex Calcium Looping reactor is given in figure HJ 
where the wall and mass thermal couplings are indicated. 
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Figure 1: In the Endex configuration the calciner temperature is lower than the carboniser 
temperature. Heat transfer by sorbent mass flow is assisted by direct wall heat transfer, 
and the CO2 pressure in the calciner is low. Light carbonation only is required; this and 
the lower calcination temperature sidestepps the sintering problem. 



Thermal coupling of reactive systems creates additional nonlinearities in 
the dynamics. In the Endex reactor configuration for CaO Looping, CO2 
separation can be considered as a gas switch in which the CO2 pressure 
in the calciner, for example, can be used to control the CO2 output mass 
flow. However, reactor coupling can amplify perturbations with adverse ef- 
fects, and a more detailed analysis of the Endex configuration is required that 
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deals with its response to perturbations, and the start-up and shut down pro- 
cesses. This analysis is performed in this paper. In section [5] we review the 
properties of Endex systems, and outline the rationale for stability analysis of 
thermoreactive systems generally and of the proposed Endex configuration in 
particular. The methodology used is described in section [3j and the dynam- 
ical model is presented. In section 0] the results are presented and discussed 
for: standalone carboniser, Endex carboniser/calciner with sorbent cycling 
but without wall thermal coupling, Endex carboniser-calciner with both sor- 
bent cycling and wall thermal coupling, and transient analysis (start-up and 
shutdown dynamics). We summarize these results in section |5j 

2. Carbonation, calcination, and Endex principles 

The chemisorption of carbon dioxide onto calcium oxide is highly exother- 
mic, releasing 170kJ/mol CO2 at 1073 K. Housed in a reactor where ther- 
mostatting is necessarily imperfect such a reaction may become thermally 
unstable and either self-quench or self-heat uncontrollably. Both situations 
are obviously undesirable from safety and economic considerations. 

Traditionally, reactors housing exothermic reactions are fitted with ther- 
mostatted cooling jackets, but the problem of maintaining thermal stability 
becomes much more difficult (and expensive) as the reactor is scaled up be- 
cause the heat removal rate is linear in the temperature and scales as L 2 
while the reactive heat generation rate is exponential in the temperature and 
scales as L 3 , where L is a characteristic reactor dimension. This scale- up 
problem is particularly relevant to reactive flue gas CO2 capture, because 
the emitter is typically a large fossil-fueled power station. 

Reaction heat that is removed by the cooling system in traditional chem- 
ical reactor plants is either dissipated into the environment, or sometimes, in 
large and sophisticated plants, partially recovered by an indirect route and 
used for another purpose. 

It is important to appreciate the potential for adiabatic, or thermally in- 
sulated, operation of an Endex configured reactor. A conventional chemical 
reactor is almost never run adiabatically (unless it is a bomb calorimeter) 
because for most common industrial reactions the adiabatic temperature rise 
for full conversion is dangerously high. Purpose-built cooling systems are 
usually necessary, which can become technically very elaborate, and expen- 
sive, when the reactor is large. The concept of an ideal adiabatic Endex 
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reactor, where the additional conservation condition of enthalpy flux conser- 
vation holds to a good approximation, suggests the very appealing possibility 
that scaling problems may be eliminated, in the same way that surface-to- 
volume ratios are irrelevant in the single conventional adiabatic reactor, while 
good conversion, thermal safety, and direct recovery of reaction heat are also 
achieved. 

An Endex system is a coupled nonlinear dynamical system, and therefore 
has the potential to exhibit thermal instabilities and is capable of more com- 
plex behavior than a single insulated or thermostatted exothermic reaction 
system. However, in the works of Gray and Ball [l|, 0] it was shown that a 
general Endex-configured reaction system can operate stably, autothermally, 
and economically, achieving almost full recovery of chemical energy and co- 
production of valuable products. The tradeoff is that the stability regime 
must be mapped for specific Endex systems in order to avoid or control (or 
even, perhaps, exploit) the coupled relaxation oscillator dynamics that are 
inherent to coupled dynamical systems. 

In the proposed Endex reactor the exothermic carbonation reaction 
CaO + CO2 <^ CaCO"3 is thermally coupled and kinetically matched with the 
reverse reaction, the endothermic calcination reaction 
CaC03 <— ^ CaO + CO2. To show that, in principle, this system can be 
operated effectively within a large margin of thermal safety we have carried 
out a linear stability analysis on the steady-state solutions of a dynamical 
model for the reactor. 

Stability analysis is a valuable reactor system design tool, and an essential 
step in the design of thermoreactive systems. Typically we are interested in 
running a continuous-flow thermoreactive system at a particular steady state, 
or set-point. Stability analysis tells us whether small perturbations around 
the steady state — which are inevitable in any real system — will decay and 
settle back onto the steady state, or grow in amplitude leading to thermal 
oscillations or uncontrollable thermal runaway. When stability analysis is 
carried out over a range of the design and operational parameters of the 
system it is often called bifurcation analysis, and it can provide a valuable 
stability map of the system. 

The mathematical theory and methodology behind stability analysis are 
well-established and accepted. Briefly, at each steady state solution of the 
parent dynamical system we construct a dynamical system for the perturba- 
tion to that steady state. Since the perturbation is small we can write it as 
a Taylor series expansion and retain only the first-order terms. According to 
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the theory of linear differential equations, the solution can be written as a 
superposition of terms of the form e Xj t where {Xj} is the set of eigenvalues of 
the Jacobian matrix of coefficients of the linearized perturbation system. A 
nonzero complex part of an eigenvalue contributes an oscillatory component 
to the solution. If the real part of an eigenvalue is positive the perturbation 
must grow exponentially with time. A stable steady state, therefore, is one 
for which all the eigenvalues of the Jacobian of the linearized perturbation 
system have negative real components. 

3. Methodology and dynamical model 

3.1. Methodology 

The following general procedure for stability analysis of thermoreactive 
systems was applied: 

1. Write down dynamical coupled mass and enthalpy balances for the 
reaction system. 

2. Choose a primary bifurcation parameter, usually a parameter that can 
be tuned experimentally, such as an inlet temperature or flow rate. 

3. Compute the steady state solutions as a function of the bifurcation 
parameter. 

4. Evaluate the stability of each solution by computing the characteristic 
eigenvalues of the linear perturbation at each point. 

5. Flag each change in sign of the real parts of the eigenvalues. These are 
the singular points. 

6. Compute the amplitude, period, and stability of any periodic solutions 
as a function of the bifurcation parameter. 

7. If appropriate, repeat the analysis using another bifurcation parameter. 
If singular points are found compute their loci using a second parameter 
to obtain a stability map of the dynamical system. 

This modelling and analysis was carried out in the well-stirred approxima- 
tion. Results from this analysis provide essential guidance for design and 
operation of an economically and safety optimised reactor system, and may 
be used to inform expensive convective simulations that require substantial 
high-performance computational resources. 
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3.2. Endex- coupled carboniser-calciner: dynamical model 

Enthalpy summation and mass balances for the gas-phase reactant lead to 
the following dynamical equations for the Endex carboniser-calciner system: 

V&i-jf = V 1 (-AH)v l (T 1 , Pl ) + FA, g (r 1)in - T x ) + (F S C S + L cx ) (T 2 - T x ) 

(2) 

dc 

V 2 -^- =V 2 v 2 (T 2 ,p 2 ) - F 2 c 2 (3) 

W 2 -^ = V 2 AHv 2 (T 2 ,p 2 ) - F 2 C 2tg T 2 + (F S C S + L ex ) (7\ - T 2 ) . (4) 

The quantities and notation are defined in table EH Appendix A. Equations 
flU)-® describe an Endex system that is heat-coupled through cycling of 
loaded and unloaded sorbent between the carboniser (denoted by subscript 1) 
and calciner (denoted by subscript 2), and via direct common wall transfer. 
The overall system is treated as fully insulated. In reality it is expected that 
heat losses to the environment will be small, less than around 5%. 

The solid-gas surface reaction rates are functions of temperature Tj, pres- 
sure pi, and fractional surface coverage 9i, (i = 1,2), 

Vl (T 1 ,p 1 ) = ^-l)6 1 ek(T 1 )C 1 S (5) 
v 2 (T 2 ,p 2 ) = (l - jg-) (1 - 6 2 ) ek(T 2 ) C 2 S, (6) 

where the rate constant fc(Tj) has the usual Arrhenius temperature depen- 
dence, and &i is given by the Langmuir isotherm 



Pi 

Pi.cq 



i+ 



Pi 



1/2 

(7) 



1/2 • 



In this model, the isotherm is based on the C0 2 molecule occupying two 
surface sites and the saturation pressure of C0 2 is pinned to the equilibrium 
pressure of C0 2 in the sorbent. 
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3.3. Data source 



Physicochemical data and reaction rates for the surface reaction system 
CaO( s ) + C02( g ) ^ CaC03( s ) were taken from the published literature and 
tabulated in |20j . The reactor design parameters and solids and gas flow rates 



are those given in [20| for a demonstration unit for scrubbing the emissions 
from a 5MW lignite fuelled power plant. Numerical values of data and 
quantities used in this analysis are given in table IA.1[ Appendix A. 



4. Results and discussion 

Standalone carboniser 

The first step in the analysis of the Endex dynamical system, equations 
CO)-©, is to analyse the carboniser segment in standalone mode. The dy- 
namics of this system will be a subset of the Endex dynamics. The standalone 
carboniser is described by the following smaller, simpler dynamical system: 

-V 1 vi{T 1 ,p l )+F 1 {c l)ia -c 1 ) © 

Vii-AIfjvtiT^px) + F 1 C ltg {T 1 ^ - Ti) + F S C S (T sM - T x ) . 

(8) 

Steady state solutions of equations ([T]) and (JSJ), using equation [51 were com- 
puted as a function of the gas residence time, T\ = V\J F^ and are rendered 
in figure |2] in terms of the dynamical variables pi (a), C\ (b), and T\ (c). The 
computations were carried out for four values of the solids flow rate F s as in- 
dicated. The real parts of the two eigenvalues of the linear perturbation were 
positive at each point, thus the steady states are stable within this regime. 

In (a) the equilibrium partial pressure, pi, eq , corresponding to each F s 
has also been plotted as dotted lines. We see that the condition p\ ^> p\ m 
holds over the range of residence times considered, thus carbonation proceeds 
spontaneously. 

In (b) we see that higher mass flow rates of solid sorbent allow improved 
uptake of the CO2 by the sorbent at shorter gas residence times. For exam- 
ple, to achieve a CO2 concentration of ~7mol/m 3 (corresponding to a CO2 
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Figure 2: Steady state analysis of the standalone carboniser. Ti t - m = 1060 K, T s ,- m = 
1021 K. 



conversion of 71%) requires a residence time of ~7.2s at a solids flow rate of 
lOkg/s, but only 4 s if the solids flow rate is 20kg/s. 

However, if it is desired to maintain the temperature below ~1123K 
and achieve appreciable CO2 conversion we see from (c) that gas residence 
times > 6s and solids flow rates > lOkg/s are required. 
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4-2. Endex- coupled carboniser-calciner: analysis 

In the standalone carboniser the heat carried by the partially carbonated 
sorbent was discarded into the environment. In effect the standalone car- 
boniser loses the heat generated by the adsorption reaction to a heat bath 
held at the constant temperature T Sj i n at the rate indicated by the last term 
in equation (jSJ). 

In the Endex-coupled system, modelled by equations ([I])-®, the heat 
generated by adsorption and carried by the partially carbonated sorbent 
is recovered directly to drive the calcination of the sorbent. There is no 
external thermostat or heat bath; instead, the dynamical temperature of the 
carboniser T\ and that of the calciner T2 are coupled at the rate indicated 
by the last term in equations <^ and 

Zero wall heat exchange 

Assuming in the first instance that L cx = 0, i.e., the carboniser and 
calciner communicate thermally only via transfer of the sorbent, we have the 
CO2 residence time in the calciner r 2 = V2/F2 as the only additional tunable 
parameter. 

The first task is to select an optimum range for the gas inlet temperature, 
T l in . The putative set-point for T l in is around 1023 K. In figures E] and @] 
steady state solutions of equations[[Hl]are plotted with T 1>in as the bifurcation 
parameter. The eigenvalue analysis gave the steady states as stable over this 
range. 

In figure |3] the calciner gas residence time T2 igas is set at 30 s, and and in 
figure H] T 2ig as is set at 60 s. The steady state solutions have been computed 
and plotted for two values of the carboniser gas residence time Ti jgaB . 

It is interesting to observe that the longer carboniser gas residence time 
Ti,gas (blue lines) gives much improved CO2 uptake by the sorbent over the 
shorter T lgas (green lines), yet the dynamical carboniser temperature Ti is 
depressed. One normally expects that an exothermic reacting system will be- 
come hotter for higher conversion of the reactant. The occurrence of the back 
reaction (calcination) at higher conversion of the CO2 should work against 
this, but the positive enthalpy change associated with the back reaction is 
not built into equation HI In this case it is the high rate of thermal transport 
provided by the mass flow of sorbent with high thermal capacitance that al- 
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lows this "temperature inversion" . When the sorbent mass flow rate F s is set 
to zero, the carboniser and calciner are decoupled and the normal increase 
in T\ with longer gas residence time occurs. 




Figure 3: Left subfigures: carboniser, right subfigures: calciner. T2 jgas = 30 s. 
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Figure 4: Left subfigures: carboniser, right subfigures: calciner. T2 !gas = 60s. 



What this means is that the sorbent mass flow is an important control 
tool for the system. We now know, for example, that T l gas can be made 
as long as you like without risking overheating of the carboniser, provided 
F s is maintained above some critical rate. In figure |5] the partial pressure 
steady states are compared for two values of F s . The higher sorbent flow 
rate evidently gives better performance. 

Comparing figures [3] and H] we note another interesting effect of thermal 
coupling of the reactors: the carboniser "sees" the calciner gas residence time 
r i,gas- Since the negative enthalpy change of the back reaction (carbonation) 
is not built into equation HJ it is the longer Ti )gas that results in warmer sorbent 
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entering the carboniser at T 2y which in turn increases T\. This improved 
thermal coupling decreases T1 — T2, consistent with the rigorous ideal Endex 
result jlj]: lim^o \T\ — T 2 \ = 0, where a = 1/F S . 

For both figures, in the carboniser segment the condition p\ 3> pi )Cq holds 
and in the calciner segment the condition p 2 <C p 2 , e q holds. The solutions 
are stable over the computed regime. However, dynamical stability — i.e., the 
behaviour of complex conjugate pairs of eigenvalues, or the occurrence of 
Hopf bifurcations — is governed by the volumetric specific heats. A system 
with low thermal capacitance is likely to become thermally unstable, and in 
fact Hopf bifurcations do occur in this system if artificially low specific heats 
are used. 

The steady state solutions with n as bifurcation parameter are plotted in 
figure O for four solids flow rates. In the carboniser a residence time greater 
than 10 s is required to achieve a CO2 uptake of 90% or more. However, 
residence times longer than about 15 s do not lead to appreciably more uptake 
ofC0 2 . 

Non-zero wall heat transfer 

In this scenario the carboniser and calciner share a common wall or walls, 
which are also heat exchangers. For cylindrical reactor segments the car- 
boniser cylinder would be embedded in the calciner segment. (The volume 
displacement is negligible and has not been corrected for in this analysis. 

For the purpose of simulating the effects of wall heat exchange we now 
regard the combined heat transfer coefficient L ex as a tunable bifurcation 
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Figure 6: Left subfigures: carboniser, right subfigures: calciner. Endex steady states as a 
function of carboniser gas residence time Ti )gas - r 2,gas = 30 s, L cx = 0. 



parameter, and the steady state solutions are plotted in figures and El The 
theoretical limit as 1/L ex — > is T\ = T2. However a reasonable upper limit 
of L ex for this configuration would be around lOkW/K. It will be shown in 
subsection C below that nonzero wall heat transfer is desirable from safety 
considerations. 

Start-up and shut-down dynamics 

Start-up: The trajectories plotted in figure M indicate that full steady 
state operation can be achieved in less than 60s, whiler quasi steady-state 
operation of the carboniser segment is achieved in less than 5 s. The time 
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Figure 7: Blue lines: T lgas = 10 s; Magenta lines: T lgas = 15 s. T\ is plotted with solid 
lines, T 2 with dashed lines. T 2:gas = 30 s, F s = 20kg/s. 
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Figure 8: (a) Ti, gas = 15 s, (b) ri !gas = 20 s. T2, gas = 30 s, F s — 40kg/s. T\ is plotted with 
blue lines, T 2 with magenta lines. 



lag is due to the much higher effective activation energy for the calcination, 
and it could be reduced by employing operational strategies that achieve 
improved effective kinetic matching of the reactions, for example, increasing 
the pumpout rate F 2 . 

Shutdown: In a normal shutdown scenario the carboniser-calciner shut- 
down dynamics are coupled to the shutdown dynamics of the flue gas gener- 
ator. Assuming the C0 2 partial pressure p c in at the inlet is reduced linearly 
and quasistatically we may simulate this type of shutdown as a quasi steady- 
state procedure, using p c in as the bifurcation parameter. In figure [TD] the 
carboniser partial pressure p\ ebbs slowly until p c ^ n is about 0.025 MPa then 
drops off dramatically. The calciner gas pressure P2 declines smoothly as the 
loading of the sorbent declines. 
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Figure 9: Cyan coloured lines: T2 jga s = 10 s. Indigo coloured lines: T2 igas = 10 s. In both 
cases Ti igas = 15 s, F s = 20kg/s, L cx = 0. In each case the initial CO2 concentrations are 
zero and the initial temperatures are set equal to the steady state temperatures. 
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Possible danger: For an Endex carboniser-calciner reactor system with- 
out wall thermal contact the only channel for heat transfer to the calciner 
is via the sorbent. If the sorbent flow F s is interrrupted (perhaps by me- 
chanical failure) the cooling capacity of the calciner is unavailable. The 
only channel for heat removal from the carboniser is then via the outflow of 
scrubbed flue gas. The adiabatic temperature rise for complete conversion of 
24.28 mol/m~ 3 of CO2 to carbonate is around 650 K. Although a tempera- 
ture rise of this magnitude would not in practice occur, because the scrubbed 
gas outflow would continue to remove heat, the scenario in which F s is in- 
terrrupted is of concern. It is simulated in figure [TTJ in which steady state 
operation is interrupted by switching off the solids flow. The time series 
were computed for zero wall heat transfer, L ex = 0, and for L ex = 1, 5 and 
10 kW/K. The thermal excursion is large and dramatic for L ex = 0, more 
than 80 K and possibly exceeding the safety limits of the vessel. Where the 
calciner and carboniser have efficient wall thermal contact, however, the ther- 
mal excursion is much smaller and quite manageable. In the case of L ex = 
10 kW/K the maximum temperature increase of 15 K occurs 100 s after the 
solids flow is switched off, after which the temperature declines slowly. 

Thus the inclusion of efficient wall thermal contact in the Endex reactor 
design may be an important safety consideration. 
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Figure 11: At t = 100 s the solids flow rate F s is switched from 40kg/s to zero, n, 
15 s, r 2 , gas = 15 s, Ti. in = 1060 K. 
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D Instabilities in the system 

In the subsections above we found that within a broad range of the ex- 
pected normal operationg conditions the Endex carboniser-calciner reactor 
is free of instabilities. This raises the question of where, exactly, in the pa- 
rameter space may instability occur in this system, since there are no global 
constraints on the nonlinear dynamical system, equations ([TJ-fllJ) that forbid 
the occurrence of a positive real eigenvalue part. 

In fact multiplicity and hysteresis may occur at low temperatures of the 
gas inlet. In figure [£2] the steady states have been computed and plotted 
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473 673 873 

gas inlet temperature, K 

Figure 12: Green curves are the stable steady states, red curves are the unstable steady 
states. F s = 5kg/s, n tSas = 2.4 s, T 2 , gas = 15 s, L cx = 0. 

as a function of the gas inlet temperature T l in . The eigenvalues were also 
monitored, and the stability of the steady states is colour-coded. 

One may, in principle, carry out an experiment in which the gas inlet 
temperature is increased quasistatically from, say, 473 K. On this lower green 
curve the reaction temperature T% increases slowly but small perturbations 
decay. When the gas inlet temperature reaches the value marked A on the 
figure the solutions change character. A real eigenvalue passes through zero 
and becomes positive. Small perturbations now grow in amplitude. A rapid 
transition to the high-temperature stable branch occurs, indicated by an ar- 
row on the figure. The inlet temperature may then be tuned up to its normal 
set-point at 1060 K, or it may be tuned down along the upper green curve. 
When it reaches the point marked B on the figure an abrupt temperature 
collapse must occur, since the solutions again lose stability. 
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Although this hysteresis loop may be exploited in start-up and shut-down 
procedures, it can be seen that it is not relevant over the normal operating 
regime of the reactor. 

Thermal oscillations associated with nonzero complex components of the 
eigenvalues do not occur in this system. This is because such behaviour 
is typically governed by the dynamics of the exothermic reaction, i.e., the 
left hand side of equation (j2J). The high thermal capacitance provided by 
the solids fraction effectively damps any ocillatory components, or complex 
parts of the eigenvalues. 

5. Summary and conclusions 

A dynamical system model was derived for an Endex coupled carboniser- 
calciner housing the CaO/CaC0 3 surface reactions, in the well-stirred fully 
insulated approximation. 

Steady state and stability analysis of the carboniser compartment in stan- 
dalone mode provided a subset of states from which to begin the more difficult 
task of analysing the full Endex system. Approximate lower bounds for the 
gas residence time T\ and solids flow rate F s were obtained. In the steady 
state and stability analysis of the full Endex system it was observed that the 
sorbent mass flow rate is an important control parameter, because the ther- 
mal transport from carboniser to calciner provided by the sorbent depresses 
the carboniser temperature at long gas residence times. 

The equilibrium partial pressures in carbonsiser and calciner segments 
were monitored. Over the parameter regime studied px jeq < p\ and p2, oq > P2, 
as required by the surface reaction thermokinetics for spontaneous carbona- 
tion and calcination respectively. 

The parameter regime studied is given in the following table: 



T-i 



l,in 



973- 1273 K 
0.1-20s 



s 



10-40kg/s 
15-60s 



T-2 



ex 



0-100kW/K 
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The real parts of the eigenvalues of the solutions remain negative over the 
parameter regime studied, hence steady state operation is stable over these 
ranges; i.e., perturbations to the system set-point decay rather than grow. 

Start-up time to steady state operation was found to be less than 60 s 
for initial CO2 partial pressure of zero and initial temperatures equal to the 
steady state temperatures. Quasi steady state operation of the carboniser 
segment can be achieved in less than 5 s. The shut-down dynamics of the 
reactor were modelled as the gradual decline in pressure of the gas inflow. 
The carboniser partial pressure fell slowly until the temperature became too 
low to sustain the reaction, whereupon the partial pressure drops off rapidly. 

Interruption of the sorbent flow was flagged as a possible source of dynam- 
ical thermal instability, in a reactor with no or poor wall thermal coupling. 
It was found that a reactor design with significant wall heat exchange be- 
tween the carboniser and calciner compartments insured against temperature 
surges in the carboniser in the event of an interruption to the solids flow. 

To complete the stability analysis, the location of instabilities in this En- 
dex system was pinpointed as a region of thermal multiplicity over tempera- 
tures well below the normal operating regime. It was noted that oscillatory 
instabilities are fully damped by the high thermal capacitance provided by 
the sorbent. 

These modelling and analysis results confirm that the proposed reactor 
configuration for the Calcium Looping reactor is a stable Endex configura- 
tion that can, in principle, scrub CO2 from a flue gas stream efficiently and 
regenerate a pure stream of CO2 for geosequestration without additional en- 
ergy requirements. The system has the potential to capture more than 90% 
of the C0 2 from flue gas emissions and release it in a pure stream, in a ther- 
mally safe reactor that requires no thermostatting. The system is stable to 
perturbations and exhibits gas switching during start-up. 
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Appendix A. Notation 



Table A.l: For compactness the subscript i is used as appropriate, 
where i = 1,2; 1 refers to the carboniser and 2 refers to the calciner. 



Ci 


concentration of CO2 


mol/m 3 


Cl.in 


inflow concentration of CO2 


Pcin/RT hin = 24.3 mol m" 


m) 


rate constant 


114exp(-i;/i?Ti)mol/(m 3 ! 


Po 




4.147e06MPa 


Pc.in 


inlet partial pressure CO2 


MPa 


Pi 


partial pressure CO2 


MPa 




equilibrium partial pressure of CO2 


p e^>(-\AH\/RT,) MPa 




reaction rate 


molrn^s^ 1 


c 1 


weighted volumetric specific heat of carboniser con- 
tents 


160 kJ K^m- 3 


Cl,g 


weighted volumetric specific heat of gas 


5.8 kJ K _1 m- 3 


C_2 


weighted volumetric specific heat of calciner contents 


25 kJ KrW 3 


^2,g 


volumetric specific heat of calciner gas 


25 J K^m" 3 


C s 


specific heat of sorbent 


975^-^-* 


E 


activation energy for calcination 


205 kJ/mol 
Vi/ti nrV 1 


Fi 


volumetric flow rate of gas into the carboniser 


F 2 


volumetric outflow rate of gas from the calciner 


V 2 /t 2 m 3 s~ 1 


F s 


mass flow rate of sorbent 


kg s _1 


L C x 


heat exchange rate coefficient 


kWK- 1 


R 


gas constant 


8.314 J/ (mol K) 
5e07m 2 /m 3 


S 


surface area 


Ti 


reactor segment temperature 


K 


7\.in 


temperature of gas at the inlet 


K 




internal volume of carboniser 


7T x 0.25 2 x 12 = 2.356 m 3 


v 2 


internal volume of calciner 


7T x 2 2 x 12 = 150.8 m 3 


AH 


reaction enthalpy 


-170kJ/mol C0 2 


e 


porosity of nascent lime 


0.51 


n 


gas residence time 


s 


Oi 


fractional sorbent surface coverage 




Ci 


carboniser solid fraction 


0.5 


C2 


calciner solid fraction 


0.008 
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